The estimate of kinetic powers of jets in FRII radio galaxies: 
existence of invisible components? 
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ABSTRACT 



We investigate the total kinetic powers (Lj) and ages (tage) of powerful jets in 



in 

Q-i- FR II radio galaxies by comparison of the dynamical model of expanding cocoons 

O ' with observations. We select four FR II radio sources (Cygnus A, 3C 223, 3C 

'^ ■ 284, and 3C 219), for which the mass-density profiles of intracluster medium 

^ ■ (ICM) are known in the literature. It is found that large fractions > 0.02 — 0.7 of 
the Eddington luminosity (LEdd) are carried away as a kinetic power of jet. The 

> ■ upper limit of estimated 2Lj/LEdd are larger than unity (< 10) for some sources. 



suggesting a possibility of super-Eddington mass accretions. As a consequence 
0\ • of the large powers, we also find that the total energy stored in the cocoon (Ec) 

exceeds the energy derived from the minimum energy condition for the energy of 
(^ ■ radiating non-thermal electrons and magnetic fields (E^i^): 4 < Ec/Emm < 310. 

This implies that most of the energy in cocoon is carried by invisible components 

such as thermal leptons (electron and positron) and/or protons. 



j_| ' Subject headings: radiation mechanisms: non-thermal — X-rays: galaxies 



radio continuum: galaxies — galaxies: individual (Cygnus A, 3C 223, 3C 284, 
3C 219) 
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INTRODUCTION 



Relativistic jets in active galactic nuclei (AGN) are a fundamental aspect of plasma 
accretion onto supermassive black holes (SMBHs). Although the formation mechanism of 
relativistic jets remains a longstanding problem, it is well established that they carry away 
some fractions of the available accretion power in the form of coUimated beam (e.g., Begelman 
et al. 1984 for review). The total kinetic powers of AGN jets Lj is one of the most basic 
physical quantiti es characterizing the jet. A lot of authors have investigated Lj in var ious 
ways so far (e.g., iRawlings fc Saunders! Il99ll : ICelotti fc Fabianlll993l : IWillot et al.lll999l ). It 
is, however, difficult to estimate Lj, since most of the observed emissions from AGN jets are 
of non-thermal electron origin and it is hard to detect the electromagnetic signals from the 
thermal and/or proton components. Hence, the free parameter describing the amount of the 
invisible plasma components always lurks in the estimates of Lj based on the non-thermal 
emissions. 

The estimate of Lj for low-power Fanaroff- Riley I (FRI) radio galaxies has been mo- 
tivated by the observations of "X-ray cavity" which is the region embedd ed in ICM with 



the s uppressed X-ray surface brightness and coincides with the radio lobe (JBohringer et al. 



19931 ). The cavities (or cocoons) are supposed to be a direct evidence of the displacement of 
the ambient ICM by the shocked jet matter. Dynamical models of cavities are a good tool, 
since the invisible plasma components as well as non-thermal electrons play a role for ex- 
pansions of cavities. For FR I sources, the total kinetic energy of the jet has been estimated 
as 



Ljtage ~ 



7c 



7c 



-PcK, 



where tage, Pc, K, 7c, a re the source age, the pressure, the volume, the adiabatic index of 
the cavity, respectively (JFabian et al.ll2002l : lAUen et al.ll2006l ). In these studies, however, the 
thermal pressure of surrounding ICM (Picm) is substituted for the one in the cavity (i.e., 
Pc ~ Licm)- This assumption may be applied only to subsonic expansions. 

On the other hand, the cocoon pressure of powerful Fanaroff- Riley class II (FRII) radio 
galaxies is expected to be larger than that of the surrounding ICM, which is expressed as 
L*c > Licm (Begelman & Cioffi 1989 hereafter BC89), and the cocoon of FR II radio galaxies 
is likely to be expanding super-sonically. Then the substitution of ICM pressure for the 
cocoon pressure is not justified. A new estimate of Lj f or FR II radio galaxies by use of the 
dynamical model of cocoon expansions is proposed by iKino fc KawakatuI (120051 ) (hereafter 
KK05), in which Lj and tage are derived from the comparison of the cocoon model with the 
morphology of the cocoon obtained by radio observations. It should be stressed that P^ is 
not assumed but solved in this model. Hence it can be applied even to the cocoons with 
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Pc > -PicM- So far, however, this estimate of Lj has been apphed only to Cygnus A. The 
expansion of the number of samples is evidently crucially important for exploring general 
characteristics of the powerful AGN jets. For this purpose, we apply the method of KK05 
to other bright FR II radio galaxies, for which the physical conditions of the associated ICM 
have been estimated in the literature. 

In the present work, we especially focus on the ratio of Lj/L^dd, where LEdd is the 
Eddington luminosity of AGN, since Lj/L-^dd is a more fundamental quantity than Lj from 
the point of view of the jet formation physics. Another interesting quantity we examine in 
this work is the ratio of the internal energy deposited in the cocoon [Ec) to the minimum 
energy (-Emin) obtained by the minimum energy condition for radiating non-thermal electrons 
and magnetic fields (e.g., Miley 1980). Some of the previous works studying this quantity 
(e.g., Hardcastle & Worrall 2000; Leahy & Gizani 2001) reported that the cocoon pressure 
expected from the inferred -Emin is smaller than the pressure of ambient matter, suggesting 
the difference between Ec and -Emin- due to the lack of minimum pressure against the pressure 
of ambient medium. Although these studies obtained the lower limit of the ratio, it is the 
value of Ec that is of greater importance. 

According to a large sample of galaxies collected recently, the frac tion of AGNs in all the 



galaxies is suggest ed to be ~ 20 — 40%, larger than previously thought (JKauffman et al.ll2003 



Miller et al.ll2003l ). The interest in AGNs is gaining momentum in the context of the co- 
evolution of galaxies and their central black holes (e.g., Kawakatu et al. 2003, Granato et al. 
2004; Di Matteo et al. 2005). AGN outflows in particular are likely to be a key ingredient in 
this context (Silk & Rees 1998; Fabian 1999; King 2003). The AGN feedback by the outflows 
may be also promising to explain the tight correlations between the ratio of the mass of SMBH 
(Mbh) to that of galactic bulge (Kormendy & Richstone 1995; Magorrian et al. 1998) and 



the ratio of Mrh and the stellar velocity-dispersion in the bulge (JFerrarese fc MerrittI 12000 



Tremaine et al.ll2002l ). In this sense, a robust estimate of the basic quantities such as L^/L^dd 



and tage of radio loud AGNs at low z is an important first step for understanding the AGN 
feedback processes in the universe. 

The outline of the paper is as follows. In §2, the model of the expanding cocoon by 
KK05 is briefly reviewed. In §3, we explain how to extract the key quantities from the 
observations of four nearby FR II radio galaxies, Cygnus A, 3C 223, 3C 284, and 3C 219, 
which are required for the comparison with the model. We then estimate the total kinetic 
power, Lj, and the dynamical ages, tage, in §4. Finally in §5, we summarize our results and 
discuss some implications on the physics of AGN jet. Th roughout the paper , we adopt a 



cosmology with Hq = 71 km s ^, flu = 0.3, and ^a = 0.7 (ISpergel et al.ll2003l ). 
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2. COCOON MODEL 

2.1. Basic equations 

Based on BC89 and KK05, we briefly summarize the cocoon model we employ in the 
following. We focus on the cocoon expansion in the over-pressured regime, namely Pc > Pa, 
where Pc and P^ are the pressures of cocoon and ambient ICM, respectively. We approxi- 
mately describe the expansion of cocoon by the following three equations: ([T]) the equation 
of the motion along the jet axis, ([2]) the equation for the sideways expansion, and ([3]) the 
energy equation. They are expressed, respectively, as 

^ = p.ir^)vlit)A^it), (1) 

Pc{t) = Pair,) Vc{t)\ (2) 

dEJt) , dVAt) 

-^ + PAt)^ = 2i. (3) 

where fj, p^.-, fh, v^,, and A^ are the velocity of jet, the density of ambient medium, the advance 
velocity of cocoon head, the velocity of sideways expansion, and the cross sectional area of 
cocoon head, respectively. Here E^, = -PcK/(7c — 1) is the total internal energy deposited 
in the cocoon, where 7c is the specific heat ratio of the plasma inside cocoon. The cocoon 
shape is approximated as a spheroid, and its volume is given as Vdt) = (47r/3)rc(t)^rh(t). 
The distance from the jet apex to the hot spot and the radius of cocoon body are obtained 
from rh(t) = [, Vh(t')dt' and rJt) = f, Vc(t')dt', respectively, and tmm is the initial 
time of source evolution. Throughout this paper, we assume 7r = 4/3, since the cocoon is 



expected to be dominated by relativistic particles (JKino et al.l 120071 ). In these equations we 
also assume that the jet has a relativistic velocity (fj ~ c) and that Lj is constant in time. 
The mass density of ICM, p^, is assumed to be given by Pa{r) = Pa{r/ro)~°', where ro and 
Pa are the reference position and the ICM mass density at ro, respectively. We set ro to be 
rh(tage), where tage is the present age of cocoon, throughout this paper. A cartoon of the 
cocoon model is illustrated in Fig. [H In this paper, we have slightly improved the model of 
BC89 and KK05 as follows: (i) a more accurate definition of Vc is employed, and (ii) the PdV 
work, which is done by the cocoon against the contact discontinuity between the cocoon and 
the shocked ambient medium, is taken into account. These corrections are necessary in the 
following quantitative estimate of Ly In fact, the estimated power is reduced by a factor of 
~ 50 from the value in KK05 for Cygnus A after taking account of the corrections (see §2.41 
for details). 
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C) -* C) 



The model parameters are Lj and t, and the unknown physical quantities are fh, v, 
and y4h. Since the number of unknown quantities is four, while that of basic equations is 
three, an additional condition is required for the system of equations to be closed. Here we 
assume that the cross sectional area of cocoon body Ac 



Tir: 



is given by A^it) oc t and 
treat X as a free parameter which is determined by the imposed condition. Once the value 
of X is determined, we obtain fh, fc, -Pc, and Ah as a function of Lj and t. It is worth noting 
that the model is capable of producing various dynamics by tuning the value of X. For 
example, the results of 2D relativistic hydrodynamical simulations by Scheck et al. (2002) 
were reproduced fairly well in Kawakatu & Kino (2006). 
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Fig. 1. — A cartoon of the employed model. The relativistic jet from FR II radio galaxy 
interacts with ICM with a declining density. Most of the kinetic energy of jet is deposited 
in the cocoon, which is then inflated by the internal energy. 



2.2. Analytical solution 



We assume that the physical quantities have a power-law time-dependence of the form 
A = A (t/tage)^, where Y is the power-law index. Then the index and coefficient of Vc, for 
example, are determined as 



Vc{t) = ^^c T 



t 



-age 



0.5X-1 



A'J' 



"age 



c (^ 

''age 



0.5X-1 



(4) 
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From this relation and Eqs. ([I])-©, we obtain the following expressions: 



-a / , ^ X(l-0.5a)-2 



^'0 



-age 



(5) 



VhiLj,t) 



^^aJ^l f — 



p^v^Ac \vqJ \t, 



X(-2+0.5q)+2 



(6) 



MLi,t) 



— _9 '-'ah 



u 



X(Q-2)(-2+0.5a)+3o-4 



(7) 



where C^h = 0.75(7^ - l)(0.5X)-"[3 - (2 - 0.5a)X]/[X{-l + 0.5a){% - 1) + 87^ - 2], 
C,c = 0.5X/7ri/2, Cpc = (0.5X)", and C^h = [X(-2 + 0.5a) + 3]-", and Vo = rh(4ge)/tage 
corresponds to the head speed assumed to be constant in time. The difference in C^h obtained 
here from that in KK05 is due to the correction made in Eq. ([3]). We assume the conditions 
of 0.5X > and X(— 2 + 0.5a) + 3 > 0, which ensure that the contribution at tmm to the 
integrations of rh and r^ is small enough. The cases that we focus on in §3 clearly satisfy 
these conditions. 



2.3. Determination of X 

As mentioned in §2.11 an additional condition that determines the free parameter X 
is required for the system of equations to be closed. In the pioneering study of BC89, 
A}^(t) = const was assumed. However, as c an be confirmed from numerical simulations (e.g., 
Cioffi fc BlondinI Il992l : IScheck et al.l |2002| ) , it is obvious that this condition is unlikely to 
hold for long-term evolutions from pc to Mpc scales. In this paper, we consider the following 
two conditions, which seem more reasonable and equally possible, in determining X: 



Constant aspect ratio (Case I) 

The aspect ratio of cocoon, TZ = Tc/ty^, is assumed t o be constant in time. This corre- 
sponds to the widely-discussed self-similar evolution (JBegelmanlllQQa : iKomissarov &: Falk 
19971 : iKaiser fc Alexanderlll997l : iBicknell et al.lll997l ). Since the time dependence of TZ 
is given by TZ{t) oc t[^(2-5-o.5Q)-3]^ ^g obtain X = 6/(5 — a) in this case. 

Constant opening angle (Case II) 

The opening angle, 9 = tan~^(y4j^ /Tr^^^rh) (see Fig. [1]), is assumed to be constant in 

time. Although there has been no previous work that has employed this condition, it 
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seems to be reasonable when the jet is precessing with a constant pitch angle. Since the 
time dependence of tan^ is given by tan^(t) oc t[0-25X(a-4)2+i.5a-5]^ X = (20 - 6a)/(4 - 
a)^ is obtained for this case. 

It should be emphasized here that these two independent conditions lead to the solutions 
that describe quite similar dynamical evolutions as long as the range of a listed in Table 
[1] is adopted (1 < a < 2). This can be seen as follows. If the constant opening angle 
is imposed (Case II), the evolution of the aspect ratio is obtained as 'R-if) oc t(2-")/(4-a) . 
It is easy to confirm its very weak time dependence, since the power-law index is limited 
to the range < (2 — «)/(4 — a)^ < 1/9. Just in the same way, if the constant aspect 
ratio is adopted (Case I), we find again a very weak time dependence of the opening angle, 
tan6'(t) oc t("-2)/[2(5-a)]^ ^i^j^ ^j^g power-law index being -1/8 < (a - 2)/[2(5 - a)] < 0. 

As a consequence of this small difference between the two solutions, the corresponding 
values of X also show only a slight difference. For example, when a typical value a = 1.5 is 
taken, we obtain X = 12/7 ~ 1.71 and X = 1.76 for Case I and Case II, respectively. Hence, 
the estimated Lj and tage based on these two sets of solutions also do not vary much, either. 
Indeed, for the given values of r^, A^, and TZ, the derived Lj and tage only differ by a factor 
of ~ 1.7 and ~ 0.83, respectively, between the two cases for a = 1.5. It is worth noting that 
for a = 2, which corresponds to the to 3C 219 (Tabled]), the two conditions give the same 
value of X = 2 and, as a result, the solutions are identical. Since only a slight change is 
found between the two cases, we focus on the widely-discussed self-similar solution (Case I) 
in the following. 



2.4. Improvements from KK05 

In the present study, we have improved the energy equation given in KK05 for more 
accurate quantitative estimations of Lj and tage- As mentioned in §2.H the resultant change 
in the derived Lj turns out to be rather large. Here we explain the reasons for this discrepancy 
more in detail. 

As mentioned in §2.H we have (i) modified dVc/dt and (ii) included the PdV term 
in Eq. ([3]). As for (i), the main flaw in KK05 is the fact that they did not take into 
account the sideways expansion in the growth of K. In fact, they employed the equation, 
dVc/dt = 27rr^t>h, whereas a more accurate expression is dVc/dt = (4/3)[7rr^t;h + 27rrcrhfc], 
which is adopted in the present study. As for (ii), it was assumed that all injected energy 
is converted to internal energy (namely, Ec = 2Ljtago) in KK05. it is obvious, however, 
that a part of the injected energy is consumed for expansions and the PdV work should be 



included, particularly for quantitative estimations. As found in §2.21 these corrections are 
reflected only in the numerical factor Cvh in Eq- © . The value of Cvh is reduced by a factor 
of ~ 3.5 owing to (i) and another factor of ~ 2 due to (ii) and, hence, by a factor of ~ 7 as a 
whole. For a given geometry of cocoon {r^, A^, and 71) and ambient density profile (pa and 
a), the derived power and age scale with the numerical factor as Lj oc C^^ and tage oc C~^, 
respectively. As a result, KK05 overestimated Lj by a factor of ~ 50. On the other hand, 
tage was underestimated by a factor of ~ 0.14, which led to the overestimation of E^ by a 
factor of ~ 14, since the latter is obtained as E^ = 2Ljtage in KK05. It is also worthy to 
note that in the present study the following relation holds: Ec — L]ta.gc- The difference of 
the factor ~ 2 arises from the fact that about a half of the injected energy is used for the 
PdV work. 



3. EXTRACTION OF THE KEY QUANTITIES FROM THE 

OBSERVATIONS 

In determining Lj and tage, we essentially follow the same procedure taken in KK05. In 
this section, we explain in detail how to extract the key quantities utilized in the model from 
the observations. 



3.1. ICM quantities 

As for the m ass-density profiles (pa) an d pressures (P a ) of IC M, we adopt t he values given 



i n the literature (JReynolds k. FabianI (1996):ISmith et alj (l2002f) for Cygnus A. ICroston et al. 



fl2004h for 3C 223 and 3C 284. and lHardcastle fc Worralll JiQQqI ) for 3C 219). In these papers, 
the X-ray surface brightness distribution of ICM w as fitted by the isothermal /9-mode l, which 
takes the form of p(r) = Pcore[l + {f /fcorcfY'^^^'^ ( ICavaliere fc Fusco-Femiandll978l ). where 
Pcore and Tcore are the core radius and density of the ICM, respectively. Since we employ 
the density profile of pi!,{r) = Pa(r/rh)~" in our model, a power-law approximation of the 
/3-model is necessary. In the present study, we determine pa = Pa(^h) and a from the /3- 
model as follows. The determination of pa is done simply by equating it with the density in 
the /3-model at the corresponding radius r^, namely pa = Pcore[l + (^h/''^core)^]~^^^^- In the 
case of Th ^ Tcore, it is clear that a can be approximated by 3(3. Only Cygnus A satisfies 
this condition, though. For the rest of the sources, rh is comparable to rcore: ^core ~ 340 
kpc, 210 kpc, and 90 kpc for 3C 223, 3C 284, and 3C 219, respectively. It is obvious that 
the above approximation of a ~ 3/3 would cause an overestimation of density gradient for 
these cases. Instead a should be taken to be a typical value in the ICM region swept by the 
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expanding cocoon. Here we determine a by requiring that Psir) should coincide with the 
density in the /3-niodel at r = 0.5rh in addition to r = r\^. Although there is no compelling 
reason for the choice of r = 0.5rh, the estimations of Lj and tage are affected little by this 
uncertainty ( §4.ip . Once pa is given, P.^ is obtained by the equation of state, which is written 
as Paij^ = -^—^Pa{r), where Ta and p = 0.6 are the temperature and mean molecular weight 
of ICM, respectively, and mn is the mass of hydrogen. We adopt the temperature used in 
the /3-model, ignoring the radial dependence of Ta as usual. In Table [H we list the values of 
Pa and Pa at r = rh and a for each source. 



3.2. rh and Ai, 



In Fig. [3l we show the VLA i mages of Cvg nus A (IPerlev et al. 



1984 ). 3C 223 flLeahv fc Perlev 



19911 ). 3C 284 flLeahv et al.lll986[ ). and 3C 219 fIClarke et al.ll 19921 ) in logarithmic scale. Con- 
tours in linear scale are also displayed to determine the position of hot spot accurately. The 
overlaid straight lines that cross each other at right angle on the hot spot are the lines we 
use to measure rh and A^. For simplicity, we neglect the projection effect of rh, which woul d 
be at most a factor of a few if we believe the unified model of AGN ( lUrry fc Padovanilll995l ). 
Ah is measured as a cross-sectional area of the radio lobe at the position of the hot spot, rh 
and Ah for each source are summarized in Table [H 

It should be noted that the plasma just around the hot spot is very fresh in the sense 
that a significant synchrotron cooling is absent. Hence, the effect of radiative cooling does 
not introduce large ambiguity in the estimation of Ah. The adiabatic cooling, on the other 
hand, is not expected to cause any problem in the estimation for the following reason. Since 
the sound crossing time in the head region of the cocoon, ~ A-^ /cg, where Cg = c/\/3 is 
the sound speed, is much shorter than the age of the cocoon, we can regard the head region 
to be uniform to the lowest order. Hence, the adiabatic cooling, if any, would decrease the 
surface brightness gradually as the distance from the hot spot increases. Contrary to this, 
the observed radio images show a sharp decline of the surface brightness at the outer edge, 
which is most naturally interpreted as the existence of the periphery of cocoon head there. 

Next we address the issue in determining Ah that arises from multiple hot spots. Double 
hot spots are actually observed in the radio lobes of Cygnus A (see, e.g., Carilli fc B arthe l 
1996). In determining Ah, we adopt the "disconnected-jet" model by I Cox et al.l (1l99ll ) 
to Cygnus A. The double hot spots are referred to as primary and secondary as follows. 
The primary hot spot is more compact and located in the inner part of the lobe, whereas 
the secondary spot is more diffuse and brighter and located in the outer part of the lobe. 
According to the disconnected-jet model, double hot spots are produced by the sudden 
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change of the jet-orientation, or the disconnection, which leads to the termination of energy 
supply to the original shock and the generation of a new jet. While the primary hot spot is 
produced by the terminal shock in the new jet, the secondary hot spot remains as a relic in 
the original jet. The schematic picture of the model is illustrated in Fig. [21 Since the primary 
hot spot is predicted to be much younger than the source age, we employ the position of 
the secondary hot spot to determine A]^. We will discuss this point more in detail in the 
following. 

The age of Cygnus A is roughly estimated to be tage ~ ''^h//9hsC ~ 2.0 x 10'^(/5hs/10"^)"^ 
yr, where jS^sC is the advance velocity of the hot spot. On the other hand, since the primary 
hot spot is observed simultaneously with the secondary hot spot, its age should be younger 
than the duration, t^ur, in which the secondary spot is bright, tdur is expressed as a sum of 
the time up to the shut-off of the energy supply from the disconnected-jet to the spot, tdis, 
and the cooling time, tcooi: 

^dur = ^dis + i^cool ~ niax(tdis, tcoolj- (o) 

The coohng time is evaluated as tcooi = niin(tsyn,^ad), where tgyn and tad are the synchrotron 
cooling timescale and the adiabatic expansion timescale, respectively. A typical value of tgyn 
at the hot spot is estimated as tsyn ~ 1-0 x lO^{B/lO'^G)-^/^{iy/l GHz)-^/^ yr, whilst tad is 
given by tad ~ ^hs/cs = 5.6 x 10^(rhs/lkpc) yr, where r^s is the size of the hot spot. These 
estimates lead to tcooi = tad for the secondary hot spot. On the other hand, tdis is given 
by tdis = Tdis/V], where rdis is the distance from the tail to the hot spot in the disconnected 
jet. Although we do not know rdis from observations, we can at least put the upper limit as 
'"h > ''^dis- We then obtain tdis < 2.0 x 10^(rh/60 kpc)(fj/c)~^ yr. Hence, from Eq. ([8]), we see 
that tdur is in the range ~ 5 x 10'^ — 2 x 10^ yr. From these estimates, it is obtained that the 
age of the primary hot spot only makes up a small fraction of its whole lifetime. Therefore, 
we adopt the secondary spot for the determination of Ah, which reflects the whole evolution 
of the cocoon (Fig. [3]). 



3.3. 7^ 

In contrast to rh and Ah, it is difficult to measure the aspect ratio of cocoon, TZ, from the 
VLA radio images, since the cocoon emission from the region far away from the hot spot is 
very dim at GHz frequency because of the synchrotron cooling (see Fig. [3]). It is well known 
that this cooling effect can be used to infer the age of radio galaxy (the spectral ageing 
method, e.g., Carilli et al. 1991). As one utilizes lower frequencies, however, the cocoon 
image is expected to be thicker, since lower-energy electrons have longer synchrotron cooling 
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Fig. 2. — A cartoon around the double hot spots in the "disconnected-jet" model. As a result 
of altering orientation, jet becomes disconnected and forms the primary and the secondary 
hot spot which are the location of the present and the previous terminal shock, respectively. 



times (e.g., Carilli and Barthel 1996). In fact, a few authors have reported the existence of 
prolate cocoons, based on the observations at a relatively low radio frequency (610 MHz) 
band (e.g., Readhead et al. 1996). It is mentioned, however, that little attention has been 
paid to observational features concerning the structure of cocoons so far. 

On the other hand, theoretical studies of jet propagation and cocoon formation with 
multi-dimensional hydrodynamic simulations clearly s upport the existence of cocoon for 



reasonably light beams going through surrounding ICM (IScheck et al.ll2002l : lAloy et al.ll2000 



Gomez et al.lll997l : iKomissarov fc Falldll997l : iMioduszewski et al.lll997l ) . The intensity maps 
of the synchrotron emissivities obtained in these hydrodynamical simulations well reproduce 
the double lobe structures as observed (e.g.. Fig. 10 in Scheck et al. 2002). It seems 
natural, therefore, to suppose that cocoons are commonly produced, although there is a 
room for further observational investigations. In the present study, we explore a wide range 
of 0.5 < 7^ < 1 in order to take account of the large ambiguity on the shape of cocoon. 
It is worthwhile to note in this re spect that the existence o f cocoon can be confirmed for 
Cygnus A in the Chandra image (IWilson et al.l l2000l . l2006l ) and that the obtained aspect 
ratio ~ 0.5 — 0.7 lies in the range explored in this paper. 
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4. RESULTS 

4.1. Total kinetic power and dynamical age 

The resultant Lj and tage are displayed in Figs. lU El and [71 Since most of the radio 
sources show asymmetries in the pair of lobes, we analyze each lobe independently. Three 
oblique lines in these figures are the obtained Lj and tage for different a's, and their ranges 
reflect the uncertainty in TZ. The thick solid line shows the result for the parameter set listed 
in Table [H The other two lines correspond to the results obtained by varying a by ±0.5. 
From the figures, it is confirmed that the results are rather insensitive to the value a. In 
each line, the power and age depend on the aspect ratio 7^ as Lj oc 7?.^""^ and tage oc 7^^~", 
and, therefore, satisfy Lj ex t~J,. Since a does not exceed 4 in any of the four sources, a 
lower aspect ratio corresponds to a higher power with a lower age. It should be noted that 
the uncertainty in the absolute values of pa and Pa is not crucial, since Pa is used only to 
judge whether the over-pressure condition is satisfied or not, and the dependence of Lj and 
tage on Pa IS rather weak, Lj oc pa and tagc oc Pa. The line outside the shaded region must 
be discarded, since the over-pressure condition is violated. The Eddington luminosity, LEdd, 
of each source is also shown in these figures for comparison. In Table [3l we summarize the 
allowed values of Lj and tage and other relevant physical parameters of cocoon obtained for 
the parameter set listed in Table [H 



Cygnus A is one of the vastly studied nearby FR II radio galaxies. iTadhunter et al. 



(120031 ) estimated the SMBH mass of Cygnus A as 2.5 x lO^M©, based on the gas kinematics 
in the narrow-line region. Its linear size is measured as rh = 70 kpc for the western jet and 
Th = 60 kpc for the eastern jet. From the employed values of rh and Ah, the power and age 
are obtained as Lj = (0.35 — 1.1) x 10^^ ergs s~^ and tage = 30 — 53 Myr for the western 
jet and Lj = (0.4 — 2.6) x 10^^ ergs s"-^ and tage = 19 — 47 Myr for the eastern jet. No 
significant difference is seen between the two jets, and we interpret that the actual age lies 
in these ranges. Note that Lj is decreased and tage is increased from those in KK05 by the 
improvement of Vc and the inclusion of PdV work. In fact, while Lj = 1.3 x 10^^ ergs s~^ and 
^age = 2.6 Myr were obtained in KK05, we find Lj = 2.6 x 10^^ ergs s^^ and tage = 19 Myr 
in this paper when the identical values of rh = 60 kpc, A^ = 150 kpc^, and TZ = 0.5 are 
employed. 

3C 223 has radio lobes that exte nd up to rh = 340 kpc in both sides. Its SMBH 



mass is estimated as 1.4 x 10^ Mq by IWoo fc Urryl (J2002l ). based on the observed stellar 



velocity dispersions. As can be seen in Fig. [31 3C 223 has asymmetry in A^. While a well 
developed cocoon head is seen at the northern hot spot (Ah = 4300 kpc^), the cocoon head 
at the southern hot spot is quite compact {A^ = 1800 kpc^). Reflecting this asymmetry. 
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the obtained power and age show quite large difference: Lj = (0.15 — 2.9) x 10 ergs s~ 
and tagc = 140 — 610 Myr for the northern jet, and Lj = (0.71 — 2.0) x 10^^ ergs s~^ and 
tage = 330 — 560 Myr for the southern jet. 

3C 284 shows asymmetry both in rh and A^. While r-^ and Ah in the western lobe are 
estimated as 380 kpc and 6200 kpc^, respectively, the corresponding values for the eastern 
lobe are 260 kpc and 4600 kpc^. The obtained power and age are Lj = (0.1 — 3.6) x 
10^^ ergs s-i and tage = 100-630 Myr for the western jet and Lj = (0.03-1.8) x 10^^ ergs s~^ 
and tage = 32 — 260 Myr for the eastern jet. Since there is no estimate of the SMBH mass 
of 3C 284 in the literature, we derive the mass from t he B-band magni t ude o f the buldge 
estimated in lShi et al.l (120051 ). By using the equation in iMarchesini et al.l ( 20041 ) which gives 



the correlation of the B-band magnitude with the BH mass, we obtain Mbh = 8.2 x IO^Mq. 

In the case of 3C 219, we only analyze the jet on the western side, since the eastern lobe 
shows severe deformation (see Fig. [3]). We could not determine Ah on the eastern side frorn 



its morphology. The central SMBH mass for 3C 219 is estimated by IMarchesini et al.l (120041 ) 
as 6.3 X 10*^ Mq. Th and Ah of the western lobe are measured as 210 kpc and 5000 kpc^. 
From these values, the kinetic power and age are obtained as Lj = (0.26 — 4.3) x 10^^ ergs s~^ 
and tage = 37 — 150 Myr, respectively. 

Large asymmetry between the pair of lobes is observed especially in 3C 223 and 3C 
284, and 3C 219. Since it seems natural to suppose that the jet properties are intrinsically 
symmetric and the power and age are identical on both sides, we expect that the asymmetry 
in the pair of lobes is due to the asymmetry and/or inhomogeneity in the ICM density 
profiles. Although this is an interesting subject, a further pursuit is beyond the scope of the 
present study. Here we assume that the actual values of Lj and tage are lying in the ranges 
obtained from both lobes. 

In Table [3] (column 4), the total kinetic powers of jet normalized by the corresponding 
Eddington luminosity, 2Lj/LEdd, are displayed. It can be seen that 2Lj/LEdd takes quite 
high values ranging from ~ 0.02 to ~ 10. In exploring the physical relations between the 
accretion power and the outflow, the total kinetic power of jet normalized by the Eddington 
luminosity, 2Lj/LEdd; is one of the most fundamental parameters. We will return to this 
issue in 55.1. 



4.2. Total internal energy vs. minimum energy 

It is intriguing to compare the internal energy, E^, deposited in the cocoon wi th the 



widely discussed energy, L^mim obtained from the minimum energy condition (e.g., iMiley 
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1980l ). Ec is linearly proportional to the total energy injected in the cocoon (2Ljtago) and 



is approximately given as Ec — -Z^j^age- The dependence on the aspect ratio is given by 
Ec ex I/jtage oc Ivj oc T?."""^ ( §4.ip . Hence, a smaller aspect ratio (or, equivalently, a larger 
power) corresponds to a larger internal energy. -Emin is the minimum value of the total energy 
(the sum of the energies in radiating non-thermal electrons and magnetic fields) required for 
a given synchrotron luminosity and is evaluated as 



E„ 



7 



247r 



-K 



3/7 



R 



12vrVV(«R)(^^fn'"^ - ^^L-"^)^"^^. 



-|4/7 



ergs, 



(9) 



where Vr is the volume of the emitting region, or is the spectral index of the synchrotron 
emission, and L^, is the synchrotron luminosity measured at frequency u, and i/min and z/max 
are the lower and higher cut-offs in the synchrotron emission, respectively, /(or) is a function 
of spectral index or which is given as 



/("r) 



3.16 X 10^2(0.145)°«(2aR + l)r(2i + 1) 



(2aR-i)r(^ + ^)r(2i + i)r(^ + |)' 

where F is the Gamma function (see, e.g., Longair 1994). 

The values of the spectral index or at the low freq uency band (178- 7 50 M Hz) and the 
flux density, E^, at 178MHz are taken from Table 1 in iHardcastle et al.l (119981 ). From the 
employed values of F^, the synchrotron luminosities are calculated as L^, = iTidlE^, where c^l 
is the luminosity distance. Although L^ is the sum of the luminosity from lobes, jets and hot 
spots, no significant overestimate of L^, is expecte d because the rad i o emission is dorninate d 
by the lobe-component for most of FRII sources (JBridle et al.l 1 19941 : IHardcastle et al.lll998l ). 
The employed values of or, L^, and other relevant quantities are summarized in Table [2J 
Here we neglect the second term in Eq. IQ in deriving -Emin, since or > 0.5 is satisfied in all 
sources. The lower cut-off frequency z/min is taken as 10^ Hz. We will comment on this value 
in the next paragraph. As noted in §3.31 although GHz radio images do not show a cocoon- 
shape clearly and only a pair of lobes can be seen, it is known that radio irnages at lower 



frequ encies reflect the cocoon shape more closely (e.g., iReadhead et al.lll996l : ICarilli et al. 



199ll ) because of the absence of efficient radiative coolings. Since we utilize a relatively 
low frequency (178MHz) band, the volume of the emission region can be approximated as 
Vr ~ K. Here we employ the median value of TZ, i.e. TZ = 0.75, in evaluating V^. We define 
the ratio of Ec to -Emin as rjc: 



Ec 



Vc 



E„ 



(10) 



The obtained values of £^min, E^ and rjc are summarized in Table [31 We find that Ec is larger 
than -Emin and rjc is in the range of 4 < r/c < 310. This implies that there is a substantial 
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deviation from the minimuni energy condition. We will discuss this topic more in detail in 



Lastly, we comment on the reliability of rj^. The lower cut-off frequency, //„ 
of the ingredients, which introduce uncertainties in rj^ because it is difficult to determine 
by radio observations. The value employed above is obtained from the following relation 
t'min ~ 10'* (i?/10~^ G)(7c,min/10)^ Hz, whcrc 7e,min is the minimum Lorentz factor of non- 
thermal electrons. Note that the the resultant -Emin does not change significantly by the 
uncertainty in i/min because of its weak dependence. For example, when a typical value 
of Or = 0.8 is employed -Emin oc (z/min/10'^ Hz)^/^^. It should be also mentioned that the 
difference between the actual emission volume and the employed one, which we do not expect 
to vary by orders, does not affect our result, since the dependence of -Emin on the emission 

3/7 

volume Vr, is weak, -Emm oc V^ . Hence the precise determination of the latter is not 
necessary. 



4.3. On the estimation of Lj and Ec 

4-3.1. Upper limits and lower limits 

It is important to consider the validity of the over-pressure condition (i.e. Pc > Pa), 
since the lower limits of Lj and Ec are determined by this condition in most cases (see Figs. 
IHlTj). In our model, a larger P^ corresponds to a smaller TZ. Though we explore a wide 
range of 7?. (0.5 ~ 1), it is intuitively more likely that TZ is smaller than unity, that is, the 
cocoon is prolate rather than spherical. Hence, the results of our analysis suggest that the 
sources examined in the present study are likely to be over-pressured indeed. Incidentally, 
the prolate shape of cocoon is endorsed b y the fact that independ ent age estimations of 



Cygnus A (ICarilli et al.lll99ll ) and 3C 284 (I Alexander fc Leahyl 119871 ) based on the spectral 



ageing method are more consistent with the results for the aspect ratio of 0.5 than for 1.0. 
(Unfortunately, the age estimations are not available for 3C 223 and 3C 219 in the literature). 

In all cases, the maximum values of Lj and E^ correspond to the minimum value of 7^ 
(Lj oc TZ'^"'~^). Though 7^ = 0.5 is chosen as the lower limit in the present study 



the possibility of even smaller aspect ratios cannot be ruled out, since the radius of cocoon 
body Tc cannot be constrained very well from the radio images. It is emphasized again that 
smaller values of 7?. predict larger Lj and Ec. 
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4-3.2. Kinetic energy in the cocoon 

We have so far neglected the kinetic energy of bulk flows in the cocoon, assuming that 
the internal energy is dominated over the kinetic energy. This may be justified by the fact 
that the cocoon is filled with shocked jet matter, which are expected to flow subsonically. 
It is, however, important to estimate the possible changes in Lj and E^ that the inclusion 
of the kinetic energy in Eq. ([3]) (energy equation) may make, since our results are rather 
sensitive to the changes in the energy equation (see the discussion in §2.4p . Although the 
lack of our knowledge on the mass deposited in the cocoon makes it difficult to estimate the 
kinetic energy quantitatively, the changes in Lj and E^ are not significant for the conclusion 
of the paper even in the case, where the kinetic energy is comparable to the internal energy, 
as shown shortly. 

Just as in §2.41 the modification in energy equation is refiected in the value of the 
numerical factor Cvh of Eq. ([6]). Note that Cvh depends linearly on the ratio of the total 
internal energy to the total injected energy, e = Ec/{2Ljt^ge), which is ~ 1/2 in the present 
study, since roughly a half of the injected energy is consumed for the PdV work. The 
inclusion of the kinetic energy modifies e as e ~ 1/(2 + /), where / is defined as the ratio of 
the kinetic energy to the internal energy and Cvh is reduced from the value obtained in the 
present study by a factor of ~ 2/(2 + /). It is hence found that even if the kinetic energy 
is as large as the internal energy, namely / ~ 1, Lj(oc C^^) and Ec{(x C^h) are reduced only 
by factors of ~ 4/9 and ~ 2/3, respectively. The inclusion of the kinetic energy, therefore, 
does not change the conclusion of this paper. 



5. SUMMARY AND DISCUSSIONS 

In this paper we have investigated the total kinetic power and age of powerful PR II 
jets. We have selected four PR II radio galaxies (Cygnus A, 3C 223, 3C 284, and 3C 219), 
for which the surrounding ICM densities and pressures are known in the literature. Below 
we summarize our main results. 

(I) Large fractions of the Eddington power in the range of> 0.02 — 0.7 are carried away 
as a kinetic power of jet in the FR II sources. 

(II) The energy deposited in the cocoon, Ec, exceeds the minimum energy, -Emin, by a 
factor O/4-310. 

Although our results allow a wide range of Lj and E^, interesting implications can still 
be obtained and will be discussed below. In §5.1., we address some issues concerning the 
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ratio of Lj to LEdd by referring to the studies of X-ray binaries. In §5.2., the energetics in 
the cocoon is constrained from the obtained rjc- 



5.1. Lj/L-Edd 

Postulating that the relativistic jet is powered by the release of gravitational energy 
(Lace) of accreting matter (e.g., Marscher et al. 2002), the jet power can be expressed as 
2Lj = ejLacc, where ej is the efficiency of energy conversion (0 < Cj < 1). Hence, 2Lj/LEdd 
gives the minimum mass accretion rate normalized by the Eddington mass accretion rate. 
Our results then suggest that quite high mass accretion rates, at least above 0.02LEdd, are 
required to produce FRII radio sources. Moreover, since 10 > 2Lj/LEdd ^ 0.65 is predicted 
for 3C 219, some FRII radio sources may have super-Eddington mass accretion rates. The 
theory of accretion disk predicts that the accretion disk of these objects is optically-thick and 
called the slim- disk (e.g., Abramowicz et al. 1988). From three distinctive X-ray properties 
(the large photon index F > 2, rapid variability and soft X-ray excess), narrow line Seyfert 
1 galaxies (NLSls) are considered to be super-Eddington objects and, therefore, are inferred 
to have a slim-disk (e.g.. Pounds et al. 1995; Boiler et al. 1996; Mineshige et al. 2000; 
Collin & Kawaguchi 2004; Shemmer et al. 2006). If 3C 219 is indeed super-Eddington, it is 
expected to show the above mentioned X-ray features like NLSls. Note, however, that 3C 
219 has a relatively hard X-ray spectrum with F = 1.58 < 2 (Shi et al. 2005), and thus the 
physical state of the accretion disk in 3C 219 could be different from those in NLSls. 

In order to explore the nature of the accretion disk in 3C 219, we compare the char- 
acteristics of AGNs with those of X-ray binaries (XRB), since both of them have common 
physical processes such as disk accretions, relativistic jets, and quenching of these jets (e.g., 
Heinz & Sunyaev 2003; Ho 2005; McHardy et al. 2006). Thanks to their much shorter dy- 
namical timescales, XRBs in various states have been observed in great detail and are found 
to occupy particular X-ray spectral states (Fender et al. 2004; Remillard & McClintock 2006 
for a review) as follows; 

(i) Lacc/^Edd < 0.01 (low/hard state: LS), 

(ii) 0.01 < Lacc/^Edd < 0.3 (high/soft state: HS), 

(iii) Lacc/^Edd > 0.3 (very high state: VHS). 

State (i) is usually accompanied by a jet. For Lacc/^Edd > 0.01, the radio emission 
is quenched in state (ii), while in state (iii), the soft VHS, which has an X-ray spectrum 
dominated by a steep power-law component (photon index F > 2), is radio quiet and the 
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hard VHS and/or the transition from the hard VHS to the soft VHS is accompanied by 
relativistic ejection events. 

If the analogy between XRBs and AGNs holds, NLSls may be in the soft VHS, since 
many NLSls have a steep power-law component and high Eddington ratios > 0.3 (e.g., Collin 
& Kawaguchi 2004), that is, higher than the upper limit set by the stability of the Shakura 
& Sunyaev (1973) disk (Shakura & Sunyaev 1973). On the other hand, 3C 219 (FR II) may 
correspond to the hard VHS and/or the transition state because it has hard X-ray spectra 
and high Eddington luminosities. In order to judge whether AGNs are scaled-up XRBs, it is 
essential to confirm that radio loud AGNs actually have the states analogous to the spectral 
states (especially VHS) in XRB^j. 

The co-evolution of a central black hole (BH) and its host galaxy together with AGN 
feedbacks have been intensively studied in various ways (e.g.. Silk & Rees 1998; Di Matteo 
et al. 2003; Granato et al. 2004). The intensive s urveys of QSOs show that the number 



density of QSOs is peaked at z ^ 2 (JFan et al.ll200ll ). The existence of QSOs at 2; > 2 with 



a central BH of a mass smaller than predicted from the bulge BH-mass relation is naturally 
expected in the build-up of SMBHs toward z ^ 2. For the exploration of the co-evolution 
processes, dusty spheroidal galaxies (e.g., galaxies emitting sub-millimeter radiations and 
ultra-luminous infrared galaxies) are the ones to be scrutinized, since the dusty-gas in them 
is one of the key quantities for the co-evolution. Kawakatu et al. (2003) pointed out the 
possibility of radio loud AGNs at high-z being QSOs in the early evolution phase (we call 
them proto-QSOs) that contain a growing BH. It is, however, difficult to explore the nature 
of proto-QSOs at high- 2; by observations in optical and X-ray bands owing to severe dust- 
absorptions. In contrast, the estimate of Lj and tage presented in this work is applicable even 
to high-z sources (e.g., Schmidt et al. 2006) without being suffered from the dust extinction. 
Therefore, the estimate of Lj/L-^dd based on the cocoon dynamics is a potential new powerful 
tool to discover proto-QSOs among high- 2; radio loud AGNs. 



5.2. The energetics 

Lastly, we discuss the energetics in the cocoon. Summing up in advance, the total 
internal energy of invisible components such as thermal leptons and/or protons tends to be 
larger than those of radiating non-thermal electrons and magnetic fields. E^ is expressed by 
components as Ec = {Uc + f/e + Uirw)Vc, where f/e, Ub, and f/jnv are the energy densities 



^It has been well established that low-luminosity AGNs are the high-mass counterpart of XRBs in LS's 
(Ho 2005). 
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of non-thermal leptons (electrons and positrons), magnetic fields, and invisible particles, 
respectively. The energy ratio rjc = E^/ E^^^ can be then expressed as rjc = {U^ + f^B + 
Uinv)/Uram, where f/min = -Emin/K IS the minimum energy density. From the obtained values 
of rjc, we investigate here the contribution of f/inv to the total energy by evaluating f/e/f/min 

and Ub/U rain- 
It is useful to express f/e/f/min and t/B/f^min in terms of Us/Ub, since Us/Ub has been 
intensively investigated by a lot of authors (e.g., Isobe et al. 2002; Kataoka et al. 2003; 
Croston et al. 2004; Kataoka & Stawarz 2005; Croston et al. 2005). Since the synchrotron 
luminosity L^, is proportional to UeU-J Vc, we obtain the relation Ue oc {Us/Ub)^^'^ , or equiv- 
alently Ub cx: (C/e/%)~'^^^ for fixed values of Lj, and Vc- From this relation and Eq. (Q, we 
can derive the following expressions: 

_^^oJ£A'^' i^^05f^V^' (11) 



Hence, f/inv/f^min is given by 




-4/7" 



(12) 



Recent observations show that the ratio of f/c to Ub is 1 < U^/Ub ^ 10 on average. Sub- 
stituting these values in Eq. flTTl) . we find 0.13 < f/B/f^min ^ 0.5 and 0.5 < Uc/Umm ^ 1-3. 
Since the obtained range of r^c is ~ 4 — 310, f/inv/f^min is evaluated as 3 < f/inv/f^min ^ 310 
from Eq. flT^ . Thus, we conclude that the internal energy of invisible particles must be 
larger than those of radiating non-thermal leptons and magnetic fields (f/g, Ub ^ f^inv) to 
explain the result obtained in this paper that Ec is larger than Emm by a factor of 4 at least. 
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Fig. 3. — Logarithmic-scaled 5-GHz VLA map of Cygnus A (upper-left) and 1.5GHz VLA 
maps of 3C 223 (upper-right), 3C 284 (lower-left), and 3C 219 (lower-right) with linearly 
spaced contours are displayed. The straight lines overlaid in each map denote the lines we 
have used to measure rh and A^. 
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Table 1: The quantities measured from observations 



Source 


Th (kpc) 


Ah (kpc2) 


Pa (g cm 3) 


Pa (dyne cm 2) 


a 


Ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Cygnus A E 


60 


150 


8.3x10-27 


8.0x10-11 


1.5 


1,2 


Cygnus A W 


70 


150 


6.6x10-27 


6.4x10-11 


1.5 


1,2 


3C 223 N 


340 


4300 


5.5x10-28 


1.2x10-12 


1.6 


3 


3C 223 S 


340 


1800 


5.5x10-28 


1.2x10-12 


1.6 


3 


3C 284 E 


260 


4600 


4.0x10-28 


6.4x10-13 


1.0 


3 


3C 284 W 


380 


6200 


2.3x10-28 


3.7x10-13 


1.4 


3 


3C 219 W 


210 


5000 


1.0x10-27 


1.6x10-12 


2.0 


4 

















Note. — Column (1) shows the names of radio sources, and the following alphabet distinguishes the pair of 
jets (see Fig. |3]). Columns (2) and (3) display, respectively, the cocoon lengths and the cross sectional areas 
of cocoon head measured from Fig. [3] . Columns (4) and (5) give the estimated ICM densities and pressures 
at r = rh. In Column(6), the estimated power-law indexes of the ICM density are presented. References for 
the density profiles and pressures are listed in column (7) 



References. 



(120041 ) :(;4 1Hardcastle fc Worralll (Il999f ) 



11 [Reynolds fc FabianI |l996h 



(2) ISmith et al. 



(200 



i; (3) 



Croston et al. 
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Table 2: The observed radio information. 



Source 


z 


dL 


Fu 


ttR Li, 






(Mpc) 


(Jy) 


(ergs s-^Hz-i) 


(1) 


(2) 


(3) 


(4) 


(5) (6) 



Cygnus A 0.0565 249 9660 0.74 

3C 223 0.1368 635 16.0 0.74 

3C 284 0.2394 1182 12.3 0.95 

3C 219 0.1744 829 44.9 0.81 



6.2x10^5 
7.7x10^3 
1.0x10^^ 
3.7x10^4 



Note. — Column (1) shows the names of radio sources. Columns (2) and (3) display, respectively, the 
redshift and the luminosity distance calculated for the cosmology with Hq = 71km s~^, Vtyi = 0.3, and 
f^A = 0.7. Columns (4) and (5) give the values of flux d ensities and spectral indexes at 178MHz, which are 
taken from the table (Table 1) of iHardcastle et al.l (jl998l ). In Column (6) the calculated luminosity densities 
at 178MHz is presented 
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Fig. 4. — The obtained ranges of power and age of Cygnus A. The three obhque hues, which 
he closely to each other are the solutions for the power-law index of the ICM density (a) 
shown by arrow; The solid line represents the solution for the estimated power-law index (see 
Table [I]). The dashed and dot-dashed lines represent the solutions for the power-law index 
increased by 0.5 and decreased by 0.5, respectively. The shaded regions show allowed ranges 
where the overpressure condition [P^ > Pa) is satisfied. Also the Eddington luminosities are 
displayed by the horizontal lines for comparison. 
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Fig. 5.— Same as Fig. H but for 3C 223. 



1048 



73 lO'*'' 
73 

03 

"3^1 o« 



1045 





^ 


= ] 


.c 


)E 


= = 






= 


^ 






























a 


= 0.f 




t>. 








T 












^ 


^ 








Kd 


d 


_ 








= 


= 


= 




1 a-1.5- 


= 


= 


= 
















^b.^^ 
























^ 
















_ 


^ 


^ 


^_ 


^^ 


^ 


_ 


^ 


^ 











— 


— 


zz 




^ 





— 


— 




























-13C284 






















E 













































1048 



73 104'' 
73 

U 
(U 

J-^IO-'B 



10 



100 

ta.e(Myr) 



lO'' 
1000 10 







= 


= 


= 


== 






= 


= 




























































a-0.9 


K.. 
















r 




'-'J^dd 

14- 1 










= 


= 


= 


E = 


^S^a-1.4- 


^ 


= 











— 


— 


— 


-N^j^ — 







— 


— 
















^v\.^ 


s 


a-1.9 














^W 














^ 




<■ 












^= 


^ 


= 


= = 




T^ 


^ 


^ 


= 
























— 








































■ 






-3C284W- 














■ 


b. 
















^S 


► 



























agev 



100 
ta.e(Myr) 



1000 



Fig. 6.— Same as Fig. H but for 3C 284. 
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Fig. 7.— Same as Fig. H but for 3C 219. 



Table 3. The obtained properties of the jet and the cocoon together with minimum energy of the radio lobe. 



Source 


k 


iage 


Mbh 


2Lj/LEdd 


E, 


-C^min 


Vc 




(10^6 ergs s-i) 


(Myr) 


(Mo) 




(10^° ergs) 


(10^° ergs) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



Cygnus A E 


0.4- 


2.6 


19 


-47 


2.5x10^(2) 


0.025 - 


-0.16 


6.2 


- 16 


1.4 


4.4 


- 11 




Cygnus A W 


0.35 


- 1.1 


30 


-53 


2.5x10^(2) 


0.021 - 


0.068 


6.1 


- 11 


1.4 


4.3 


-7.8 




3C 223 N 


0.15 


-2.9 


140 


-610 


1.4x10^(1) 


0.16- 


-3.2 


30- 


130 


0.88 


34- 


150 




3C 223 S 


0.071 


-0.2 


330 


- 560 


1.4x10^(1) 


0.078 - 


- 0.22 


12 


- 22 


0.88 


14 


-25 




3C 284 E 


0.3- 


18 


32- 


260 


8.2x10*^(3,4) 


0.053 


- 3.4 


26- 


210 


1.8 


14- 


120 




3C 284 W 


0.1 - 


3.6 


100 


- 630 


8.2x10*^(3,4) 


0.018- 


- 0.67 


21- 


130 


3.0 


7- 


43 




3C 219 W 


2.6- 


43 


37- 


150 


6.3x10^(3) 


0.65 


- 10 


130 


- 500 


1.6 


79- 


310 






























to 

00 



Note. — Column (1) shows the names of radio sources, and the following alphabet distinguishes the pair of 
jets (see Fig. |3]). Columns (2) and (3) display, respectively, the total kinetic powers and ages of the radio jets. 
In column (4) and (5), the black hole mass and the kinetic powers normalized by the corresponding Eddington 
luminosity are displayed, respectively. References for the central SMBH mass are given in parentheses. In 
Columns (6) and (7), the total energies deposited in the cocoon and the minimum energies required for the 
synchrotron emission are displayed, respectively. Column (8) gives the ratios between Ec and -Emm- 
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